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We present analytic solutions to a class of cosmological models described by a canon- 
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ical scalar field minimally coupled to gravity and experiencing self interactions through a 
hyperbolic potential. Using models and methods inspired by 2T-physics, we show how ana- 
lytic solutions can be obtained in flat /open/closed Friedmann- Robertson- Walker universes. 
Among the analytic solutions, there are many interesting geodesically complete cyclic solu- 



tions in which the universe bounces at either zero or finite sizes. When geodesic completeness 



is imposed, it restricts models and their parameters to a certain parameter subspace, in- 
cluding some quantization conditions on initial conditions in the case of zero-size bounces, 



but no conditions on initial conditions for the case of finite-size bounces. We will explain 
C3 ■ 

the theoretical origin of our model from the point of view of 2T-gravity as well as from the 
point of view of the colliding branes scenario in the context of M-theory. We will indicate 
how to associate solutions of the quantum Wheeler-deWitt equation with our classical an- 
alytic solutions, mention some physical aspects of the cyclic solutions, and outline future 
directions. 
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I. INTRODUCTION 



In this paper we will analytically express a cyclic universe using exact solutions in a 
scalar-tensor theory with a scalar field o minimally coupled to gravity. 
The full action of our theory is 

S = Jd 4 x^^R(g)-^d^d v a-V(a)Y (1) 



where the potential is 



Via) 



4 



b cosh 4 ( — % ) + c sinh 4 ( —= | 

\/6/ \VqJ 



(2) 

Here b and c are dimensionless free parameters of the potential, and n~ l is the reduced 



Planck mass k _1 = y — 2.43 x 10 18 ^-. A plot of the potential energy V (a) for various 
signs and magnitudes of b, c, consistent with stability (b + c) > 0, show that the profile of 
this potential is similar to those often used in the study of cosmology. This potential was 
chosen because we can solve the equations exactly, thus enabling us to perform the type 
of analysis presented in this paper. We assume that the general features discussed here go 
beyond the special choice of potential 2 . 

The model of Eqs. ffTFZl) was initially inspired by 2T physics as described in [6| 

and section HI Ah below. The same model fits also in the worldbrane scenario 171, as inspired 

fi n . n 

by D-branes in M-theory [8j. The ideas of a cyclic universe [2] modeled in Ref.[9[ can be 
adapted to reproduce the same potential V(a), thus describing a universe that consists 
of two 3+1 dimensional orientifolds that periodically collide with each other by oscillating 
in an extra fifth dimension. It is quite interesting that this connection emerged between 
2T-gravity and M-theory. In subsections (II Ap and (II Bp we will comment on the different 
origins that converged on this model. 



In a previous cosmological application of this model [6| V (a) was an energy density of 
the order of the grand unification scale (^gc/t) 4 - in that case, 6 1 / 4 v / 6/t _1 or c L l At \f§KT 1 were 



2 In fact, this is not the only potential for which we are able to give a full analysis with the complete 
set of analytic solutions [1]. In the near future we will present a similar discussion for the potentials 
Vi (a) = (^f) ^ e -2«WV6 + ce -4na/V6^j and ^ 2 = be 2 PK o ■/ <sfl ) where &j Cj p &re dimensionless 
real parameters. The profile of V\ (a) , with c > and b < 0, is similar to the profile of the potential used 
initially in the cyclic cosmology model in [2] . 
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of order tugut ~ 10 16 ^§^, thus leading to dimensionless values for the parameters b or c in 
the order of 10 -12 . Exact solutions have a way of finding applications in various fields. For 
physical applications of our solutions, including cyclic cosmology or other future cases, the 
value of the parameters b, c, should be chosen appropriately depending on the application. 

The complete set of analytic solutions for this model, in a homogeneous, spatially flat, 
isotropic Friedmann-Robertson- Walker (FRW) universe, were obtained in our earlier paper 



6|, and some of their perturbations were studied in j6](l0|. In the current paper we will 
emphasize a subset of these solutions that are geodesically complete and describe a universe 
smoothly evolving through big bang or big crunch singularities at which the universe shrinks 
to zero size, but then it continues to perform periodic expansions and contractions that 
describe a cyclic universe, all without violating unitarity or the null energy condition in a 
flat universe. We will also include the effect of spacial curvature for the FRW universes 
(k = 0, ±1) in our new exact solutions, and we will exhibit cyclic solutions with finite size 
bounces as well. 

Perturbations such as radiation is easily included in the exact solutions, while anisotropy 
can be discussed with analytic approximations; but those aspects, as well as the quantum 
treatment through the Wheeler-deWitt equation, which require more detailed discussions, 
will appear in a separate paper [11 ]. 

The complete set of homogeneous, isotropic classical solutions presented in {f| show that, 
the generic solutions for the field a (r) and the scale factor a (r) describe a geodesically 
incomplete geometry. The geodesic incompleteness can be exhibited in terms of conformal 
time t as defined by the line element ds 2 = a 2 (r) (— dr 2 + ds 2 ) , where ds 2 is the line element 
of the 3-dimensional space. As an illustration consider the spatially flat case ds 2 = dx ■ dx. 
The geodesic x M (r) of a massive particle in this flat geometry is described by its velocity 

dxjr) _ p ^ 

dr y/fp -f m 2 a 2 (r) ' 

where p is the conserved 3-momentum. In terms of this conformal time r, the complete set 
of solutions in jsj show that, for the generic solution, the scale factor a (r) starts out at zero 
size at some time a (ti) = and grows to maximum size a (Y2) = a max in a finite amount 
of conformal time (T2 — T\) =finite. It turns out that a max is infinite in the case of b > or 
finite in the case of b < 0. Furthermore a (r) has this same behavior in an infinite number of 
different disconnected intervals of conformal time r. Each such separate interval describes a 



universe that starts out with a big bang and expands to maximum size. Moreover there are 
other disconnected intervals in which the universe contracts from maximum size to zero size. 
Evidently such generic solutions of the Friedmann equations are geodesically incomplete. 

If expressed in terms of cosmic time t defined by the line element ds 2 = (—dt 2 + a 2 (t) ds 2 ) , 
the geodesic equation reads 

dt ~ a (t) yjp 2 + m 2 a 2 (t) 

where a (t) is expressed in terms of cosmic time t. The relation to conformal time is dt = 
a (r) dr or t (r) = J'^a (r) dr, where t (n) = defines the big bang at a (t (ti)) = 0. Hence 
t (r) is given by the area under the curve in a plot of a (r) versus r, for some interval 
t\ < t < T2, starting with the big bang. An example of a geodesically incomplete curve 
a (r) is Fig.l in [6j while examples of geodesically complete ones are given in many figures 
in the current paper. Since a (r) in the generic solution is given in disconnected r intervals, 
the cosmic time t (r) cannot be defined for negative values before the big bang. Hence the 
geodesic equation above is artificially stopped at the big bang at the finite value of time 
t (n) = 0. This is one of the signs of geodesic incompleteness of this geometry. In addition, 
when the area under the curve is finite (the solutions for b < 0), the total cosmic time 
t (r 2 ) is also finite, and geodesies are again artificially stopped at a finite value of cosmic 
time t (t 2 ) , providing another sign of geodesic incompleteness. In this way, each interval 
that is geodesically incomplete in conformal time appears again as geodesically incomplete 
in cosmic time. This type of geometry bounded by singularities, and classical solutions in 
them, occur often in General Relativity and are commonly used in its applications, as in 
our own paper Q. But in view of the geodesic incompleteness of the generic solutions of 
the Friedmann equations displayed in the conformal frame, it feels that this must be an 
incomplete story. 

We think that a more satisfactory approach, especially for cosmology, is to find those 
solutions that describe a geodesically complete geometry. This type of solution is what we 
will describe in the current paper. It turns out that among the classical solutions presented 
in 6] there are some unique solutions that are geodesically complete in the Einstein frame. 
In this solution the patches of conformal time in which a (r) is real in the Einstein frame 3 are 



All solutions, including those that are geodesically incomplete in the Einstein frame, are geodesically 
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smoothly connected from r = — oo to r = oo. Then the universe sails smoothly through sin- 
gularities, while geodesies of both massless and massive particles smoothly continue through 
singularities to the next cycle of the cyclic universe. 

The requirements for such solutions depend on whether the bounce is at zero size 
a (r ban g) = or finite size a {j bang ) > 0. In the case of zero size bounce, that occurs when the 
spacial curvature is zero, k = 0, initial conditions of the two fields need to be synchronized 
and periods of oscillation need to be relatively quantized, as we will described in detail. 
These requirements result in some quantization conditions among the available parameters 
consisting of integration constants of the differential equations for a (r) , a (r) as well as the 
parameters b, c in the potential energy V (a) , and also radiation and anisotropy parameters 
when they are included. Because of these requirements these geodesically complete zero-size 
bounce solutions are associated with a countable set of boundary conditions (but still an 
infinite set, in the sense described in section (II D[) ) . In the case of finite-size bounce none of 
these constraints occur on boundary conditions, but in this case there is spacial curvature, 
k/rQ ytz o, which needs to be large enough to compete with the potential energy. In this 
case, as long as the parameters that define the model are within a certain continuous region, 
the generic solution is the finite-size bounce solution without any further requirements on 
boundary conditions. 

When perturbations, such as radiation and anisotropy are included, or when quantum 
effects in the form of the Wheeler-deWitt equation are taken into account, there still are 
geodesically complete solutions that have a similar character to what we will present in this 
paper. They still form a countable set for zero size bounces, while they are the generic 
solutions for the finite bounces. Either way, the distinguishing character of geodesic com- 
pleteness has an appeal that seems important for physical applications, as we will discuss in 



a future paper [111 ]. 



complete in other frames. Indeed, in the 7-gauge that we will discuss in section (HU), all solutions are 
geodesically complete. However, as viewed from the point of view of the Einstein frame, a (r) for such 
solutions becomes imaginary in the patches that complete the geodesies in the 7-frame, and hence the 
physical meaning in the Einstein frame becomes obscure. We intend to study this phenomenon in a future 
paper, but for the current paper we concentrate on only the geodesically complete solutions in which a (r) 
is real for all r. 
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A. 2T-physics origin 



We would like to briefly summarize here the main points of how the model in Eq.flT]) 
relates to 2T-gravity ^ijs] in 4-space and 2-time dimensions as the conformal shadow in 
3-space and 1-time. More generally, according to 2T-physics, a theory in lT-physics in 
id — 1) + 1 dimensions is one of the many shadows that comes from d + 2 dimensions 
A useful shadow that appeals to the intuition of physicists accustomed to relativistic field 
theory is the one called the conformal shadow. In the conformal shadow there is a local scale 
symmetry (Weyl symmetry). The original 2T-gravity in 4+2 dimensions does not have a 
Weyl symmetry in 4+2 dimensions, instead this crucial gauge symmetry in 3+1 dimensions 
is dictated in the conformal shadow as a remnant of general coordinate transformations in 
the extra 1+1 dimensions 5|. Other less familiar shadows provide other descriptions of the 
physics that are related by duality transformations to the conformal shadow, and often they 
can provide hidden physical information that is systematically missed in the conventional 
formulation of lT-physics 



Besides this duality aspect, 2T-physics may also provide additional constraints on the 
interactions of fields in lT-physics even within the conformal shadow. The constraints in 
theories of interest, in the conformal shadow, are mainly on scalar fields and their interac- 
tions. These constraints have been determined generally in [14j in the presence of gravity 
or supergravity. Most of the emergent constraints can be rephrased as being consequences 
of various symmetries in lT-field theory, but not all of them. Some of those additional 
constraints are not motivated by fundamental principles in lT-physics, as discussed in jl^ . 
so they can be taken as signatures of 2T-physics. Here we will deal only with the simplest 
version of scalar fields that obey the constraints. This is the case when all scalar fields 
are conformally coupled to gravity in the IT version. This is a familiar form in IT field 
theory, but lT-gravity does not require that all the scalars be conformal scalars; by contrast 
2T-gravity has this as an outcome for the conformal shadow in one of its allowed versions 

n 

(more general form of constrained scalars in [141]). 

Moreover, in the conformal shadow there is no Einstein-Hilbert term in the action, so 
there is no dimensionful gravitational constant that would break the local Weyl symmetry 
explicitly. Instead, the Einstein-Hilbert term emerges from gauge fixing the Weyl symmetry 
within the conformal shadow (see below). This is a mechanism called "compensating fields" 
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which is familiar in conventional field theory. Such structures of 2T-physics are compatible 
with the construction of satisfactory models of a complete theory of Nature directly in 
4+2 dimensions, including the standard model 15 1, its generalizations with supersymmetry 

nn nn n 

[16(|17| or grand unification, gravity |4j[5|], supergravity [14J, all of which lead to applications 
in LHC physics and cosmology {(5]. 

The derivation of the conformal shadow in 3+1 dimensions from the 4+2 dimensional 
theory is described in detail in . For the purpose of the current paper it is possible to 
skip this detail and start out directly in 3+1 dimensions by requiring a local scale symmetry 
(Weyl symmetry). Then, the ordinary looking model in Eqs. f[TTl2"j) can be presented as a 
gauge fixed version of the following gauge invariant field theory in 3+1 dimensions, which is 
a conformal shadow that contains one scalar field s (x) in addition to a dilaton <fi (x), both 
conformally coupled to gravity as follows 



S = J d'x^-g Q^<R0 ~ \g" u d,sd u s + 1 (0 2 - s 2 ) R (g) - <ff f^j ) . 



(5) 



The field 4>(x) has the wrong sign in the kinetic term, so it is a ghost (negative norm 4 ). 
This sign of the kinetic term is required by the Weyl symmetry if the sign in front of the 
curvature term j2<f) 2 R (g) is positive. However, due to the local scale symmetry the ghost 
can be gauge fixed away, so this theory is unitary. The gauge symmetry is preserved for 
any potential of the form 4 / where / (z) is an arbitrary function of its argument 
z = |. In this action there is no Einstein-Hilbert term with a dimensionful gravitational 
constant, but instead, the factor (0 2 — s 2 ) 1 plays the role of a spacetime dependent effective 
"gravitational parameter" . 



B. Braneworld origin 

A cyclic model, inspired by D-branes in M-theory [8(, was developed in [2j where it 
was discussed for a very different potential than Eq.(T2]). However, it is possible to recover 



4 There are models of cosmology based on the notion of "quintom matter" 18| which also introduce a 
negative norm ghost field. We should emphasize that those models have actual ghosts and therefore are 
non-unitary and fundamentally flawed. Despite some similarity, our model is fundamentally different 
because of the Weyl symmetry that eliminates the ghosts, thus having fewer degrees of freedom. Our 
action, our solutions which do not violate the null energy conditions, and the discussion of the physics are 
also different. 
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precisely the current model of Eqs. ffT|2]) in the colliding world brane scenario as follows. 
One should compare Eq.(27) in Ref.[9j to the model in Eq.(j5]) before gauge fixing the Weyl 
symmetry. Both models have a Weyl symmetry that is a remnant of general coordinate 
transformations in extra dimensions (although the extra dimensions in the two cases do 
not have identical signatures). One should compare our fields here s (x) , (x) to the fields 
ip^ (x) in Ref.Q], since they are both conformally coupled scalars and have precisely the same 
kinetic energy terms. Furthermore, the potential V (a) of Eq.((2]) is also recovered, if one 
replaces the unknown terms in Eq.(27) of Ref.js], 2Wcft [g + }~ 2Wcft [<7~] + S m [g + ] + S m [g~], 
by just a cosmological term on each brane. Namely, replacing the unknown expression by, 
b + y/g+ + b-yfg~, where b± are constants, and using their definition of g ± , gives J2± b±V~g^ = 
J2± b±y/—g~ (V' ± ) 4 • This is indeed the potential 60 4 + cs 4 in Eq.([5]), which in turn leads to 



the potential V (a) after the Weyl gauge symmetry is fixed to obtain the Einstein frame as 
described in jg] and below. 



C. Fixing the Weyl symmetry 

The Weyl symmetry can be gauge fixed in several forms. In the Einstein gauge de- 
noted by a label E, such as 4>Ei $e, g 1 ^ '■> the gauge is fixed such that the curvature term 
Y2 ((p 2 — s 2 ) R (g) becomes precisely the Einstein-Hilbert term (^e) , so that in the Ein- 
stein gauge we have 

1 ^ 1 



12 V'*-*"*) = 2^ (6) 
In this gauge it is convenient to parametrize <f>E, se in terms of a single scalar field a 

0E(l) = ± ^ cosh (^M), S£W ^ sinh (^|)). (7) 

Then the gauge fixed form of the action in Eq.([5]) takes precisely the form of Eq.(pQ), where 
the potential V (a) is arbitrary as long as the function / (z) is arbitrary. 

The Friedmann-Robertson- Walker metric (FRW) in this gauge takes the form 

ds% = -dt 2 + a% (t) dsl = a 2 (r) (-dr 2 + dsj) , (8) 

ds\ = i — tttt + r2 ( dd2 + sin2 ed( i )2 ) ; k = °» ±l - ( 9 ) 

1 — kr /r 

where ds 2 is the metric of the 3-dimensional space, r is the conformal time and a (r) = 
aEt/t^)) is the cosmological scale whose dynamics we wish to study in this paper. The 



relation between ordinary co-moving time t and the conformal time is 5 



dt = a (r) dr, or t (r) = / a (r 1 ) dr'. 

'o 



The scalar curvature of the metric in Eq.flH]) is given by 

R(9e 



6 ,'d k 
+ 



a 2 a 



in 



(11) 



(12) 



where r is a constant radius that sets the scale of the curvature of 3-space 6 when the 
dimensionless scale factor is a = 1. 

Thus, in this gauge, for the purpose of homogeneous solutions of the equations of motion, 
the dynamical variables are a (r) and a (r) which interact with each other as prescribed by 
the action ([1]). Their equations of motion reduce to the Friedmann equations [19] as follows^ 



K 

y 

a 4 



2a 2 



Via) 



k 



1 9 ' 



K 
If 



O 



-Via) 



a aa Tr , . . 
- + 2— + V (a) 



0. 



(13) 
(14) 
(15) 



We had previously found all the exact solutions of these equations for the potential V (a) 
given in Eq.([2]) and a flat universe k = 0. These were tabulated in p. In this paper we will 
emphasize the subset of those solutions that are geodesically complete and in addition we 
will generalize them by including non-zero spacial curvature k = ±1. Further generalizations 
including radiation and anisotropic metrics will be given in To explain what we mean 
by a geodesically complete solution we need the following discussion. 



In this paper the overdot denotes derivative with respect to conformal time a (r) = da/dr and a (t) = 
d 2 a/ (dr) 2 . The derivative with respect to comoving time t can be rewritten by using the chain rule as 
It = o{j)1t' ^ ot example , the Hubble parameter H ~ 1 
are expressed as 



CLE {t) dt 



%f and its derivative + H 2 = -^^M , 



a(r) dH 2 

H _ 2TT ' ~JT + H 
a A (T) dt 



a 

773" 



(10) 



6 The parameter tq sets the scale for the curvature. If normalized to today's curvature, with r$ representing 
todays size of the universe, then K = k/r 2 . is extremely small even when k ^ 0. In that case we can 
completely forget the effect of spacial curvature. However, there are cosmological models that play with 
the curvature parameter as applied in the early stages of the universe. In that case ro may be within a 
few factors of 10 of the Planck scale, in which case the curvature is enormous. In order not to miss on 
possible interesting solutions we will not pre-judge the size of this term and investigate the solutions that 
emerge. Then in various physical applications we may or may not neglect the parameter K = k/r^. 
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D. Geodesic completeness 



Note that the Einstein gauge in Eq.([6]) can be chosen only in patches of spacetime 
when the gauge invariant quantity [1 — s 2 /(f) 2 (x M )] is positive. This quantity may 
be expressed in the Einstein gauge (i.e. when it is positive only) as (1 — s 2 E /(p 2 E ) = 
(cos1i(k<7/v^6)) 2 • We must note that this gauge invariant quantity could vanish at various 
times r. We did in fact find that it does vanish at various values of r in generic analytic 
solutions for a (r) ,a (r) given in [6|. However, when (1 — s 2 E /(j) 2 E ) vanishes 4>e diverges so as 
to maintain the gauge choice for the gauge dependent quantity ((f) 2 — s 2 ) as given in Eq.flHJ). 
But the Einstein gauge was fixed under the assumption that ((f) 2 — s 2 ) was positive; if it can 
vanish can it also change sign? This is the question that initially motivated two of us (J to 
study this model in the (f>, s version, rather than the a, a version. Will the dynamics require 
the gauge invariant quantity (1 — s 2 /(f) 2 ) to change sign in some patches of spacetime, thus 
creating patches with antigravity? If yes, what would that mean cosmologically for the 
universe we live in? 

In our previous study in jsj our exact solutions for (f>, s showed that generically the 
dynamics does require the gauge invariant (1 — s 2 /(f> 2 ) to change sign. However, the point 
at which (1 — s 2 / (f> 2 ) vanishes corresponds to a big bang or a big crunch singularity where 
the scale factor in the Einstein gauge vanishes a 2 (r) = a E (t (r)) = (recall a (r) is gauge 
dependent), so the physical interpretation for our own universe may be stopped exactly 
where (1 — s 2 /(f> 2 ) vanishes, and therefore the solution could be stopped artificially at that 
moment in conformal time r. This is geodesically incomplete, as well as gauge dependent 
from the point of view of a (r) as defined in the Einstein frame. But nevertheless, if one 
insists that the theory is defined only in the Einstein frame, one could accept a geodesically 
incomplete patch for a physical interpretation in the usual interpretation of gravity. This 
type of geodesically incomplete solution, which is very common in applications of gravity, 
was used in the application to an inflating universe we discussed in our previous paper 6j. 

In the current paper we take the point of view that geodesically complete solutions are 
more satisfactory. To discover and better understand the solution, we examine the factor 
(4> 2 — s 2 ) that multiplies R (g) in the action. To overcome the gauge dependent description, 
of the Einstein or other frames, we focus on the gauge invariant quantity (1 — s 2 /(f) 2 ). The 
point at which it vanishes corresponds to the big bang or big crunch. When it is positive 



11 



we can choose the Einstein gauge to describe ordinary gravity, but in patches when it is 
negative there is antigravity. Only the geodesically complete solutions has no antigravity 
by having (1 — s 2 / cf 2 ) > as a function of r. Of course, quantum corrections near singu- 
larities may smooth out the behavior of solutions. Notwithstanding the cloudiness of our 
understanding of quantum gravity at this time, it still seems to us attractive to identify the 
geodesically complete geometries and solutions in the applications to cosmology, expecting 
that this feature survives the quantum effects, as it seems to be the case at the level of 



the Wheeler-deWitt equation llfl . Hence we will identify the circumstances in which there 
are geodesically complete solutions in which the quantity (1 — s 2 /<f) 2 ) never changes sign 7 . 
We found that this is indeed possible, and we explicitly obtained those unique geodesically 
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complete solutions that are presented in the present paper and in 

As we will see in the explicit solutions given below, it turns out geodesic completeness, for 
the bounce at zero size, requires two ingredients. First, the parameters in the model have to 
be in a certain range and satisfy certain quantization conditions. In other words not every 
model can yield geodesically complete cosmological solutions with zero size bounces. For 
example, in the case of the flat FRW universe and in the absence of any perturbations, the 
ratio of the parameters b/c in the potential V (a) above must be in the range — j < (b/c) < 4 
and must be quantized as in Eq. (|30|) . This condition on b/c is relaxed in the presence of more 
parameters, such as curvature (k = ±1) or radiation, but there is always one combination 
of parameters and initial conditions that is quantized. Second, even with the quantized 
parameters, initial conditions for (r) , s (r) must be synchronized with each other in order 
to generate geodesically complete solutions in which (1 — s 2 /4> 2 ) never changes sign. If initial 
conditions are not synchronized, then (1 — s 2 / cf 2 ) will change sign and all solutions will be 
geodesically incomplete; but this is what we want to avoid, and on this basis we consider 
the solutions with synchronized boundary conditions, namely only the geodesically complete 
ones, as being those that provide a fuller story of cosmology. 

We have also found exact analytic solutions, that obey <fi 2 — s 2 > at all times, in which 
the initial conditions need not be synchronized or the parameters of the model need not be 
quantized. Such solutions occur in the presence of spacial curvature in the closed universe 
k = 1, and provide a cyclic cosmology where the universe bounces at a finite size. For this 



7 I. Bars thanks Paul Steinhart for stimulating discussions that led us to focus on this question. 
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to be possible the curvature needs to be large enough to compete with the potential energy 
V (a) , as will be discussed in sections fllVCIIIVDj) . 

Since (1 — s 2 /(f) 2 ) never changes sign for such solutions the physics at all times is com- 
patible with Einstein's gravity, since then one can indeed choose the Einstein gauge, Eq. ([()]), 
at all times in such a universe. 

If one takes the point of view that the theory is defined directly in the Einstein frame 
in terms of a, a, as in Eq.([T]), then the <p, s configurations in which (1 — s 2 /(f) 2 ) changes sign 
is a spurious outcome of the parametrization in terms of (f>, s in Eq.Q. In that case all 
field configurations in which (1 — s 2 /<f) 2 ) is negative are excluded by definition. If one also 
requires geodesic completeness then the solutions we present below are the only ones that 
satisfy the criteria. 

This avoids the question of what happens to the physics for those solutions that 
are geodesically complete in a more general sense than the Einstein frame, by allowing 
(1 — s 2 /(f) 2 ) to change sign. If the theory is defined at a more fundamental level (as in 2T- 
physics, or as in the colliding branes scenario) in which one would accept all the consequences 
of the action in the 0, s version of Eq.Q, then one must investigate the properties of those 
solutions as well. What we do know from our explicit solutions (6), in the cases in which 
initial conditions are not synchronized or parameters are not quantized, is that the quantity 
{(f) 2 — s 2 ) does not remain negative after switching sign, but oscillates back to positive. So, it 
appears that the universe recovers from antigravity and comes back to a period of time with 
ordinary gravity. However, if allowed to continue its motion in complete geodesies, the sign 
changes back and forth again and again. Perhaps the physics appears to be all wrong dur- 
ing the time periods (or more generally spacetime patches) where ((f) 2 — s 2 ) is negative, but 
we don't really know the physical cosmological consequences of such solutions for our own 
universe. We think that it would be interesting to find out eventually the physical viability 
and meaning in cosmology of the more general geodesically complete solutions allowed by 
the action Eq. (J5j) . So we will not throw away yet the generic solutions which were included 
in the list in (pj), nor will we settle the associated physics questions in this paper. So, at a 
less ambitious level, for the moment we concentrate only on the geodesically complete cases 
that also satisfy 2 — s 2 > 0, as required by the action Eq. ([!]). 
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II. ANALYTIC SOLUTIONS 



To analyze the model in the 0, s version we find it useful to gauge fix the Weyl symmetry 
in other forms. A very useful gauge is to choose the conformal factor of the metric to be 1. 
We will call this the 7-gauge. In this gauge we will denote the fields with a label 7, such 
as 7 , s 7 , g^. Then, the Robertson- Walker metric in this gauge loses the scale factor since 
a 7 = 1 

dv^ 

ds 2 = -dr 2 + -— T + r 2 (d9 2 + sin 2 Odcf 2 ) , (16) 

1 — kr /r 

and its curvature is a constant given by 

k 

R (g 7 ) = 6K, with K = — , k = 0, ±1. (17) 

r o 

In this 7-gauge there is no scale factor for the universe, but now both 7 (x) , s 7 (x) are 
dynamical variables, with 7 having the wrong sign in the kinetic term. The advantage of 
this gauge is that the dynamics of 7 , s 7 become much simpler and can be solved exactly 
in certain cases. After obtaining the solution one can transform back to the Einstein gauge 
to find a (r) , a (r) . For the case of only time dependent fields the gauge fixed form of the 
action ([3]) is 

L = \ + - \ (-^ + 4) - 0V (i) . (18) 

But one should also remember that r reparameterization symmetry of general relativity 
requires the vanishing Hamiltonian constraint (this is the Goo = ^00 Einstein equation) 

h = \ (-4 +p!) + f (-4 + 4) + f'f (J) = °. < 19 > 

where the canonical momenta are = — 7 and p s = s 7 . The negative norm ghost is 
eliminated because of this constraint on the phase space (0 7 , s 7 ,p<^,p s ). The Wheeler-deWitt 
(WdW) equation of our theory (0, s) = takes an interesting form in the 0, s basis 

Q n - di) + § (-0? + , 2 ) + 0V (i) ) ^ (0, S ) = 0. (20) 

As a side remark, note that for K > (closed universe) and in the absence of the 
potential, 4 /(|) = 0, the system in Eqs.f TTBti2T?|) describes the SO(l, 1) Lorentz symmet- 
ric relativistic harmonic oscillator in 1+1 dimensions, with (0, s) representing the ("time", 
"space" ) directions respectively. As in other cases of harmonic oscillator in several dimen- 
sions, this system has a larger hidden symmetry, which is SU(1, 1) I)SO(l, 1) in this case. 
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The quantum version of the relativistic harmonic oscillator (i.e. the WdW equation for 



/ {s/4>) 

paper 



0) was studied and solved exactly in sections VI, VII and Appendix of a recent 



20| by using unitary representations of SU(1, 1) . This may be taken as a toy model 
to begin a study of the WdW equation for our case 8 . Of course, we are interested in the full 
WdW equation, including the potential energy 4 / ( 4 ) , radiation, and anisotropy, as will 



be discussed elsewhere 



We now turn to the classical equations of motion satisfied by 7 , s 7 , including the po- 
tential energy. Such classical solutions provide a semi-classical approximation to the WdW 
equation. We are interested in an exactly solvable case so that we can study the issues we 
raised with certainty. One of those exactly solvable cases corresponds to a special form of 
the potential, namely 4 /(s/</>) = 60 4 + as 4 , that in turn corresponds to the hyperbolic 
potential V (a) given in Eq.(T5]). The equations of motion for </> 7 (r) and s 7 (r) are directly 
obtained from the Lagrangian or Hamiltonian given above. But it is also instructive to derive 
the equations directly from the Friedmann equations in Eqs.( IT3]IT5"l) by using the following 



The WdW equation Hty (<fr, s) = is satisfied by an infinite set of solutions of the relativistic harmonic 
oscillator [20(. These are \& (</>, s) = J2^Lo c nVVi (4>) V'n (s) , where the c n are arbitrary and ipn (4>) , 4>n (s) 
are the standard 1-dimensional harmonic oscillator wavefunctions that satisfy the eigenvalue equations 
\ (-d\ + K4> 2 ^ $ n ((f>) = E n if> n ((f)) and \ (~d 2 s + K s 2 ) ifj n (s) = E n ^ n (s) , where E n =VK{n+\). For 

all these solutions \& (0, s) has an overall gaussian factor exp[— ^M- ((f) 2 + s 2 )] times a polynomial. This 
shows that the probability distribution |^ ((f), s)| 2 for these generic solutions is not purely "timelike" (not 
4> 2 > s 2 ), but rather it covers both "timelike" and "spacelike" regions in (<fr, s) space. This is not surprising 
since for generic boundary conditions the classical equations also do not obey 4> 2 (t) — s 2 (t) > at all times. 
Only special boundary conditions with synchronized phases for <j> (r) , s (r) at r = can yield classical 
solutions that have this property, as we have explained in the text. Similarly, the relativistic harmonic 
oscillator has just one quantum state whose probability distribution is concentrated in the timelike region 
<f) 2 > s 2 ; has a damping factor of the form exp (— ((f) 2 — s 2 )) and vanishes on the "lightcone" <p 2 = s 2 . This 
was given in the appendix of [20j (interchange spacelike with timelike in that appendix). This solution 



is the "timelike singlet of SU(1, 1)" [20(. There is also a separate "spacelike singlet of SU(1,1)", while 
the other generic solutions correspond to a superposition of other non-singlet unitary representations of 
SU(1, 1) [20(. Referring to the comments of the last two paragraphs in section ppp . if the fundamental 
action is in the Einstein frame (as in Eq.([T])), then only the timelike SU(1, 1) singlet is acceptable as a 
solution of the WdW equation. If, on the other hand, the starting point is more general (as in action 
(J5J), then in order to favor only the solutions that are consistent with <f> 2 — s 2 > 0, the model should 
have an additional ingredient. This could be an appropriate potential energy term, effects of radiation, 
curvature, etc., or some appropriate constraint that is natural in the model. This would then characterize 
the "right" model. 
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connection between the 7-gauge and the Einstein gauge (derived in [6J]) 



K 2 

2 ™ ( ±2 „2^ 



\/61 



which gives 



*= (21) 



6^ + c 4 

Inserting these in the Friedmann equations we obtain the equations for 0, s as follows 

= 7 - 460 7 + #0 7 , (23) 
= s 7 + 4cs 7 + Ks 7 , (24) 

= - + ^) - + c, 7 + \k4) ■ (25) 



The important observation is that the second order equations for 0, s decuple from each 
other, so they are exactly solvable. The third equation simply states that the energy of the 
solution must be matched to the energy of the s solution = E s . Once the solution 
is obtained it is transformed back to the Einstein frame by using Eqs.( l2Tj) . thus providing 
the desired solutions for a (r) , a (r) in the Einstein frame. The general generic solutions of 
these equations, for all possible ranges of the parameters and boundary conditions arc listed 



in I for the K = case. The solutions are expressed in terms of Jacobi elliptic functions 9 



as given in jg] and below. 

Now we focus on the subset of solutions that satisfy the criteria we laid out. When 
the spacial curvature of the FRW universe is zero (i.e. k — 0), we found that geodesically 



9 The Jacobi elliptic functions that we need for our solutions are denoted as sn (z\m) , cn (z\m) , dn (z\m) . 
These are periodic functions that have properties similar to trigonometric functions. The following for- 



mulas 
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are directly useful to verify our solutions explicitly. The derivative of Jacobi elliptic functions 



are given in terms of expressions somewhat similar to those for trigonometric functions. 

— sn (z\m) = cn (z\m,) x dn (z\m) , (26) 

— cn (z\m) — —sn (z\ni) x dn (z\m) , (27) 

— dn (z\m) — —m x sn (z\m) x cn (z\m) . (28) 

They also satisfy quadratic relations, such as 

(sn (z\m)) 2 + (cn (z\m)) 2 = 1; m (sn (z\m)) 2 + (dn (z\m)) 2 — 1. (29) 
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complete solutions, with ((f) 2 — s 2 ) > at all times, can occur only when the ratio b/c takes 
on the following quantized values in the range — | < | < 4, with c positive, c > 0, and 



4c 



with n = 1,2,3, ••• . (30) 



Note that each value of n defines a given model. The explicit solutions are given in 
Eqs. (I32|I41||49T) for b > 0, b < and 6 = respectively. When additional parameters, such 
as curvature, radiation, etc. are included in the model, then this quantization condition on 
the model is relaxed, but still some combination of parameters and integration constants 
must be quantized as we will discuss. Furthermore, if the spacial curvature is sufficiently 
large, so that the curvature term in the action can compete with the potential, then we 
find, geodesically complete, finite, bouncing solutions of a cyclic universe where the bounce 
occurs at a minimum finite size of the universe. For such cases there is no quantization 
condition on the parameters of the model, but instead, the initial conditions on the fields 
must be within a certain range defined by those parameters. 

Before we give the mathematical details, we first explain how the physics is easily captured 
by interpreting these decupled equations in terms of an analog mechanical problem of a 
particle moving in a potential. In the case of <fr the Hamiltonian is H ((f)) = + V ((f)) , 
with V (<f)) = \K(f> 2 — b(p 4 , while in the case of s the Hamiltonian is H (s) = |s 2 + V (s) , 
with V (s) = \K(f) 2 + cs 4 . According to Eq.( l25i) the only acceptable solutions for 0, s are the 
ones that satisfy 

H((f)(T)) = H(s(r)) = E (31) 

The corresponding potentials are depicted in Figures- ( 1TTI22 112 7h for the cases of k = 0,±1. 
We have included the cases of positive b (heavy solid curve V ((f))) and negative b (dashed 
curve V ((f)))- We have drawn the pictures for only positive c (solid thin curve V (s)) while 
for negative c the V (s) curve is reflected from the horizontal axis in each figure. 

These figures, combined with the physical intuition of a particle in potential, capture the 
physical aspects of our solutions. We approach the mathematical analysis systematically for 
each figure and investigate the various ranges of the parameters b, c, K and the integration 
parameter E. We will start with the simplest case of zero curvature and analyze it thoroughly 



in section fllllj) . We will then discuss the positive/negative curvature cases separately in 
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sections (llVIIVj) respectively. 



III. THE FLAT (k = 0) FRW UNIVERSE 

We first discuss the flat case (k = 0) that has the fewest parameters. For k = the only 
possible solutions are for E s = > as shown by the horizontal dashed line in Fig.([T]). In 
the case of s, the potential is an infinite positive well (V (s) = cs A with c > 0), therefore the 
particle is trapped in the well, and s (r) oscillates back and forth between turning points 
—s (E s ) < s (r) < so (E s ) given by cSq = E s . In the case of 0, if b is negative (dashed curve 
V {(f)) = —b(f> A with b < 0) its behavior is similar to the one just described for s, so <fi (r) also 
oscillates back and forth — O (E^) < (f> (r) < O (E^) at the energy level = E s . But if b 
is positive, then the particle is in an inverted well (heavy curve V ((f)) = — b(f> 4 with b > 0). 
So, at the energy level E$ = E s , which is higher than the peak of the hill, the particle 
will come up the hill from (f> = — oo, go over the top of the hill, and slide down the hill to 
infinitely large positive values of (f>. The trip may also happen in reverse direction depending 
on initial conditions. It turns out that the trip from <f) = — oo to = oo is completed in 
a finite amount of conformal time r, so allowing all possible values of proper time, (f> (r) 
repeats the trip periodically again and again by jumping from <f) = oo to <f) = — oo. Such 
solutions (given analytically in [6]) solve all the equations but do not yet address the sign 



In order to have the oscillation amplitude of (f> to be larger than the amplitude of s it 
is necessary to have — | < b < 4c (consistent with the curves as drawn in the figure). The 



of (0 2 (r)-s 2 (r)). 




FIG. 1: The flat FRW universe, k = 0. 
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lower bound — | < b for negative b, is partially understood from the figure which shows 
that the turning point for should be at a greater distance from the origin as compared to 
the turning point for s. However the \ factor in — | < b and the upper bound b < 4c for 
positive b, emerge from the details of the solutions in Eqs.( l32|4H49l) . This is a constraint 
on the model. If the potential energy V (a) does not satisfy this property it will not be 
possible to maintain 2 (r) > s 2 (r) at all times. In addition, to insure 2 (r) > s 2 (r) , we 
must (i) synchronize the initial conditions of the 0, s particles at the origin at r = 0, namely 
(0) = s (0) = 0, and (ii) also require that their periods are commensurate, so that at the 
time r when returns back to zero s also returns to zero at the same time (although s could 
make several returns to zero in the meantime). Commensurate periods can be arranged only 
by quantizing the parameter b/c. This too is a condition on the model. If the quantization is 
not satisfied in the model (0 2 (r) — s 2 (r)) will change sign periodically as a function of time. 
But when the potential V (a) satisfies the required conditions the model yields geodesically 
complete solutions in which (0 2 (r) — s 2 (r)) never changes sign, but periodically touches 
zero, which corresponds to a big crunch smoothly followed by a big bang. This is just the 
solution we sought as given in Eqs.( |32ll4Tll49l ). We see that the model has to be "right" to 
be able to yield such a solution. 



A. b > case 

The solutions that satisfy this description are a subset of those in 6j and explicitly given 
by the following expressions. For positive b, c, the only geodesically complete solution occurs 
for the quantized values of b = 4c/n 4 , with n — 1, 2, 3, • • • , as follows 

, ( T ) - Kn sn (M) ( \ _ « sn (^\l) (32) 

" [ ) fWcTl + cn{^\\Y MTJ ^WcTdn^W) 

Here the Jacobi elliptic functions, sn (z\m) etc. (see footnote [9]) , appear only for the case of 

the parameter m = 1/2. The energy level = E s is parametrized in terms of the parameter 

T which provides a scale for conformal time, as 

E* = E.= (33) 

where T (or E s = E^) is one of the integration parameters that appears in integrating 
the differential equations. Note that this T is related to the overall factor in Eq.( )32|) that 
determines the amplitude of oscillations of s 7 (r) . 
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It is easy to verify that these are solutions of Eqs. fl23H25|) by using the properties of Jacobi 
elliptic functions given in footnote [9) The plot of these functions in Fig.(j2]) conveys their 
periodic properties and show how 7 (r) > si (r) at all times. The quantum n = 5 chosen 
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FIG. 2: 7 (r) and s 7 (t) plotted for n = 5, k = \/6, T = 1, c = 1/8. 

for this figure corresponds to the ratio of the periods of versus s. The times at which and 
s vanish together 7 (r) = s 7 (r) = are the only times when 0^ (r) = s 7 (r) , at which point 
the universe goes through smoothly from a big crunch to a big bang. At an intermediate 
time t = Ttum the quantity 7 (r) — s 7 (r) attains a maximum; this is the turnaround point 
at which the universe stops expanding and begins contracting. These features are seen in 
Fig-© for the scale factor a (r) , and the scalar a (r) which are given by Eqs. fTST]) . 




-10 - 



FIG. 3: a (r) and a (r) plotted for n = 5, k = \^6, T = 1, c = 1/8. 

A parametric plot for 7 (r) , s 7 (r) is given in Fig.(j3]), with on the horizontal and s on 
the vertical. This captures similar information to Fig.(T5]). It is for the model b/c = 4/n 4 
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with n = 5, which leads to the 5 nodes in the figure. The time after the first node is a 
fast inflation period, as seen also from Fig.(j3]). In a semi-classical approach to the Wheeler- 
deWitt equation, the curve shows the region in (0, s) space where the WdW wavefunction 
$ (<f), s) is expected to have the largest probability. This is the unique curve (for the n = 5 
model) purely in the "timelike" region <fi 2 (r) — s 2 (r) > in (0, s) space. The corresponding 
WdW wavefunction is the analog of the "timelike SU(1, 1) singlet" in footnote OH]). 



X*X 



FIG. 4: The arrow at the origin marks the crunch/bang moments and the arrows at the ends mark 
the turnarounds. 



Other quantities of interest to convey the properties of the solution include the Hub- 
ble parameter H — -% (see footnote EJ), the kinetic energy of the a field, K (r) = 
its potential energy V(cr(r)) and the equation of state parameter given by w(r) = 
(K (r) — V (t)) / (K (r) +V(r)). Their plots appear in Figs. fl5ll6ll7j) . 
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FIG. 5: Hubble parameter H (r) and a (r) for n = 5. 



The Hubble parameter decreases from infinity at the big bang, quickly approaching a 
constant at the turnaround (with a few small ripples depending on n), switches to negative 
at turnaround and then slowly reaches negative infinity at the big crunch. 

The potential and kinetic energies of the a field are fairly close to each other in magnitude 
most of the time. At the turnaround the kinetic energy vanishes K (T turn ) = while the 
potential energy is a constant V (r turn ) = 6 2 6/k 4 (since a = at turnaround). Both K (0) 
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FIG. 7: The equation of state w = (K - V) / (K + V) , for n = 5. 



and V (0) are infinite at the bang or crunch, but V (0) is larger at the singularity since 
w (0) = —1/9 as seen in Fig. ((7]). 

The behavior of various quantities near the bang/ crunch singularity is better understood 
by studying the Taylor expansion near r — > for any value of n as follows 
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The last expression shows that w = — 1/9 at the singularity for all values of n. This behavior 
seems to be surprising according to common lore. 

We emphasize that this behavior near the singularity is only for our geodesically complete 
analytic solutions that satisfy both the relative quantization of their periods as well as the 
synchronization of the initial conditions. If either of these is not satisfied (i.e. for non- 
geodesically complete solutions in only the Einstein frame) the behavior near the singularity 
is radically different. This behavior seems to be of measure zero in the space of all solutions. 



In our next paper 



11] we will further analyze this issue including the effects of anisotropy 



and the quantum effects via the Wheeler-deWitt equation. 



B. b < case 

We repeat the same type of analysis for b < which refers to Fig.([T]) with the dashed 
curve representing V (0) . Geodesically complete solutions occur only if b/c has one of the 
quantized values b/c — —1/ (n + l) 4 , with n — 1, 2, 3, • • • . Then the unique solution is 

M n + i) s " (<rfgHO , , _ « ">{t\i) (41) 
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We provide a few plots similar to the ones in the previous subsection for the model n = 5. 
Despite some similarities, there are notable differences in the behavior as compared to the 
b > case of the previous section as indicated in the following comments. 




FIG. 8: <f)(r) ,s (r) , have finite amplitudes. FIG. 9: a (r) has a finite maximum. 



From Figures 8,9 we see that the universe grows up to a maximum finite size before it 
turns around. Another way of plotting the information in Fig. 8 is the parametric plot for 
(j) (r) , s (r) in Fig.t flOl) ; note the 5 nodes corresponding to n = 5. As in the previous case, 





•A 











FIG. 10: Crunch/bang is at the origin, turnaround at the edges. 



this figure is associated with the semiclassical probability distribution of the Wheeler-deWitt 
wavefunction in ((f), s) space. 
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FIG. 11: Behavior of a (r) ,H(r) FIG. 12: Temporary inflation periods. 



From Figs. (11,12) we see that there are temporary inflation periods during which H (r) is 
temporarily almost a constant, and the acceleration 10 is positive, as seen in Fig. 12. The 
number of such temporary acceleration periods is determined by n. 

It is interesting to speculate on whether this could be a mechanism to explain the current 
accelerated inflation state of our universe; namely could it be that we currently are in such 
a period which is inflationary only temporarily on the scale of the lifetime of the universe? 




FIG. 13: Energy components of a field. FIG. 14: Equation of state w (r) . 



The acceleration is denned in the Einstein frame as d ■ This may be written in terms of conformal 
time as \d T (^d T a) . For the purpose of the plot we have defined the quantity "acc", as the acceleration 
divided by an extra factor of a. 
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The energy of the a field is small except near the bang/crunch where it is infinite. The 
equation of state w (r) grows to infinity at turnaround w (T turn ) = oo, while it takes the 
value w(0) = -1/9 at the bang/crunch where V (0) /K (0) = 5/4, while V (0) , K (0) are 
both infinite. 

The behavior of various quantities near the bang/crunch singularity is given by the Taylor 
expansion near r — > for any value of n as follows 
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(42) 
(43) 
(44) 

(45) 
(46) 
(47) 
(48) 



C. 6 = case 



Finally, for vanishing 6 = 0, the solution corresponds to the n — > oo limit of either the 
positive or negative 6 branches, and is given by 

We provide a few plots similar to the ones in the previous subsections. 
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FIG. 15: 4>(t) grows linearly with r. FIG. 16: There is a single crunch/bang. 
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FIG. 17: Parametric plot equivalent to Fig. 15. 

These plots correspond to the n = oo limit of the previous plots for either b > or b < 0. 
Therefore their interpretation is similar to the discussion above 
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FIG.18: H(t) , <t(t) decrease. 



FIG. 19: Temporary acceleration periods. 
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The temporary acceleration periods persist. 
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FiG.20: The energy of a field decreases. 



The behavior of the energy, pressure and the 
The behavior of various quantities near the 
expansion near r — > as follows. These agree 
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FIG. 21: Equation of state w (r) . 
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equation of state are indicated on the figures, 
bang/ crunch singularity is given by the Taylor 
i with the n = oo limit of the previous cases 
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IV. THE CLOSED (k = +1) FRW UNIVERSE 



For k = +1 the system of equations is Eqs. (l23H25]) . This amounts to the motion of two 
particles <ft, s satisfying the equations of motion derived from Hamiltonians 



with 



H (0) = + V b (4>) , H(s) = ^s 2 + V c (s) 



V c (s) = X -Ks 2 + cs\ and V b (0) = l -K<f - b<p\ 



as plotted in Fig. (l22|) for K = +i/tq, and whose energies are constrained by 

Ha, = H K . 



(60) 



(61) 



(62) 



The motion changes character depending on whether = E s is larger or smaller than the 
peak of the inverted double well in Fig. (122]) . i.e. the maximum of Vfe>o(0)- This critical 
value is given by 



E* 



K 2 
16b' 



(63) 



Therefore we need to discuss separately the high and low energy levels E^ = E s above and 
below this critical value as shown in Fig. ( 122]) . 




FIG. 22: The closed FRW universe, k > 0. 



A. Higher level E > E* , and b > or b < 

For the higher level of E s = Ea > E* , the intuitive physics discussion works in exactly the 



same way as the discussion for the k = case at the beginning of section (IIIII) . for both b > 
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or b < 0. In the case of b > 0, the particle s 7 (r) is trapped in an infinite well and oscillates 
between turning points ±so (E), while <p 1 (r) oscillates from minus infinity to plus infinity. 
In the case of b < both particles are trapped in infinite wells, so they oscillate between 
turning points ±</>o {E) and ±Sq (E) respectively. The turning points ±0 O (E), ±so (E) are 
the points where the curves V c (s) , Vb<o (0) intersect the horizontal curve E0 — E s — E as 
seen in the figure. 

Hence for all cases b > or b < at the higher energy level of E s = E^, the geodesically 
complete motion is described by plots of 7 (r) , s 7 (r) that are similar in character to the 
k = case given in Figs.(T2]-21). In particular for very small curvature K (large values of the 
curvature radius r ), the K ^ plots should approach the K = plots. Therefore, we will 
not include the K ^ plots here. 

The exact solutions for the higher energy level are denoted as (r) , <p+ (r) where the 
superscript "+" refers to E > E* with E s = E^ = E. The solution s+ (r) is given by 

+ ( \ l l - K2T s SU fe K ) m * ( E ) = \ C 1 - KT s (£)) ^ 
M ^ V 8cT * dn(^|m.)' T s (£) = (lQcE + K 2 )~ 1/4 1 J 

while (r) has the following expressions for b > 



1 sn | 




m+ J 


1 


867 +l + 01(^1 


m+j 



, . , , , m + (E) = \ + K 1%(E) 

-° (T) -i/8w?; r^TT^TT' T+ (B) s (6 46B)-v-. • (65) 



or 6 < 



" 7 '" V 8 l fo l T2 dn(£|m_)' T_(E)^(16|6|E + ^r 1/4 

Note that all symbols T, m that appear in the Jacobi elliptic functions are determined in 
terms of the energy level E = E s = E^ as given above. The T± (E) are determined in terms of 
E, b, c, K by the energy condition E^ = E s = E, by using the following expressions computed 
from the Hamiltonians H(<fi),H(s) given in Eqs. (I60H62]) for the solutions s+ (r) , (r) 
above 

E ° = m* (1 " K2Tt) • Et = m? E * = (1 " I<2Tt) ■ (67) 

Note also we have assumed that E s = E^ = E is higher than the critical value E* in Eq.( l63|) . 
This yields the expressions for T Sj ± (E) , m S) ± (E) given above as well as the ranges for these 
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' 1 - K 2 Tt sn (?:\ m -) m_ (E) = |(1 - KT 2 (E)) 



(66) 



parameters as a function of the energy level as follows 

KT 2 (E) < KT 2 (E) < _L f KTl (E) < ± (68) 

which determines the possible range of values for m s , m± as E changes in the range E > E* 

~ (l - (1 + c/ \b\y 1/2 ) < m s (E) <l (~ - ^) < m_ (£) < i; ~ < m+ (£) < 1. 

(69) 

It is easy to see that in the zero curvature limit K — > these solutions reduce to Eqs. (!32|) 
for b > or Eqs.flU)) for 6 < 0. 

As they stand these solutions do not yet satisfy the requirement <fi 2 (r) > s 2 (r) at all 
times. This can be satisfied only by requiring the period of to be a multiple integer of the 
period of s. The analytic expression for this conditions is 

b > : T + Q (m+) = 2nT s Q (m s ) , n = 1, 2, 3, • • • (70) 
b < : T_Q (m_) = nT s Q , n = 1, 2, 3, • • • (71) 

where the quantity Q (z) is a well known special function, namely the quarter period of the 
corresponding Jacobi elliptic function, and is given by the following integral representation 
of the EllipticK integral 

f"/ 2 ijf) 

Q(z)= - = EllipticK (z) . (72) 

</o V 1 — z sin 9 

The consequence of this is to require a certain combination of parameters (b, c, K, E) to be 
quantized. The range of parameters in the model (b, c, K, E) that can satisfy the constraint 
can be determined numerically 11 by using the expressions above. Therefore only a model 
that can satisfy this condition can give the corresponding geodesically complete solutions. 

B. Lower level E < E* , and b < 

When the energy is less than the critical value, the exact solutions are denoted as 
s~ (r) , (p~ (r) where the superscript "— " refers to the energy interval < E s = E^ = 

11 For example, using Mathematica, which recognizes the function EllipticK(z), one can plot 
T±Q (m±) /T S Q (m s ) as a function of one of the parameters E, &, c, K (example b) while the other three 
are chosen arbitrarily. When the plot matches an integer n, this fixes the value of the remaining parameter 
(i.e. b in the example above) in terms of the integer n and the other three parameters. 
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E < E*. We consider at first the case of b < 0, for which Vfe<o (0) is represented by the 
dashed curve. There is no difference in this case with the b < case above, so the formulas 
above apply, namely s~ (r) is the same as s+ (r), and 4>~ b<0 (t) is the same as 0^ b<o (t) , 
except for the fact that now the energy is in the range < E < E* 

_ \ \-K^ sn m s (E) = |(1- KT S 2 (S)) 

j i-icn™ (£1™-) m_ (£) = i (i - £:t 2 (JSO) 

07 ' b<o(Tj V ^ dn(^\m_y T_(E) = (16\b\E + Kr l/4 ' 
Since the energy is less than the critical value, < E s = E^ < E*, we must now restrict the 
range of the parameters m, T to 



< KT S (E) < 1, — < i^T_ (£?) < 1, (75) 



\b\+c~ s v 7 " ' v/2 



which implies 



0< m ,( B )<I , m _ (B) <i(l--L). (76) 



To obtain geodesically complete solutions the analog of the quantization condition in Eq. (ITT 
must be further imposed, T_Q (m_) = nTQ (m). 



C. Lower level E < E* , and 6 > 0, finite bounce 

For the case of b > (Vf, >0 (0) represented by the inverted double well in Fig. ( 122]) ) there 
are new features. For <p now there is the possibility to either be trapped inside the false 
vacuum, or be outside of it, depending on initial conditions. When <ft is trapped in the false 
vacuum it should oscillate between two turning points; when it is outside it should oscillate 
between a finite value and infinity. So the solution has the form 

/ KT 2 _ l / _ x m in (E) = [KTl (E) - 1) 

^;l>o(r)= V * sn — \m in , / / - v-V* , (77) 

V2&T in V Tin J Tm (E)=if + ^Jf- AbEj 

t y/KT** + 1 1 m out (E) = -\ (KT* ut (E) - 1) 

K ~ V2bT out m (^\ mout Y T 0Ut {E) = {K*-mET^ ' 

V out J 



32 



Meanwhile s (r) continues to oscillate as before between two turning points, so it is still 
given by the same expression, namely 



'l-K*T* m (r.\ m >) m s (E) = |(1- KT 2 (E)) 



■ s 



The energies of these solutions are computed in terms of the parameters m, T by using the 
Hamiltonians in Eqs. (16011 62]) as follows 

1 KT 2 — 1 K 2 T 4 — 1 

771 { -t r^2rp4\ jpin xl * in x jpout ±s - * out x /orA 

E ' = w^^- KT ')^ E *=-mir iE * = mr \ t ■ m 

All energies must be positive and equal to each other E^ ut = E s = E and E™ = E s = E, as 
well as smaller than E*. This yields the expressions for T Stin out (E) and m s ^ intOUt (E) given 
above as well as the ranges for these parameters as a function of E, as follows 



1 > KTl (E) > y ^qjp 2 > ^ 2 n (JB) > 1, KT 2 out (E) > 1. 
Similarly, the range of values for the parameters m s , m, n , m ou( are then as follows 



< m s (E) < X - 1 1 - y ^A-j , < m m (£7) < 1, m out (£7) < 0. (81) 

To obtain the geodesically complete solution, in the case of the inside solution the quan- 
tization condition Ti n Q (m,; n ) = nTQ (m) is required^. However, in the case of the outside 
solution no quantization condition is needed as explained below. 

We now comment on the outside solution given by 0~' 6 °>o ( r ) > s jb>o ( r ) m Eqs. ( 175179^) 
because it is different in character as compared to all the previous cases. It describes a 
periodically contracting/expanding universe with bounces that occur at minimum finite 
values of the scale factor, while the maximum is infinite. 0~' b °"o (t) describes the behavior 
of 4> outside of the false vacuum, the amplitude for <p is always larger than the amplitude 
for s^ 6>0 (r) at all times and for all boundary conditions, including the additional arbitrary 
parameter tq. This solution represents cyclic bounces at finite minimum sizes of the universe. 
This is shown in Fig. (l23|) . The cyclic bounce occurs for all values of c > 0, all values of 
b > 0, and all values of the relative initial conditions tq at r = 0. This is why we added 
an additional phase tq in the expression of 0~' fe °"o ( r )- There is no need to synchronize the 
boundary conditions at r = for (f)(0) ,s (0) in order to get geodesically complete solutions 
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FIG. 23: (A-yb^o ( r ) anc ^ s -y b>o ( r ) f° r the bouncing solution. 

that satisfy <fi 2 (r) > s 2 (r) at all times. Moreover, 6/c need not be quantized since the 
periods of <p, s can now be arbitrary relative to each other. 

The corresponding plots for the scale factor a (r) and a (r) for the bounce are given 
in Fig. (1241. The minimum size of the scale factor a (r) ~ * <\> 2 (r) — s 2 (r) at the bounce 




FIG. 24: The bounce and turnaround. 

instant is not the same each time since the initial values of (r) , s (r) are not synchronized 
and their periods are not related. 

For this bounce solution to play a physical role in cosmology we need the curvature 
terms K<ft 2 , Ks 2 to be able to compete with the potential terms 60 4 , cs 4 . For this to be 
phenomenologically tenable in a complete cosmological model, a period of inflation must 
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follow after the bang so that the universe inflates to its current size and to its almost flat 
current condition (since r$ would not be identified with today's size of the universe). 

As a limiting case of the above solutions we point out the special case of the b > scenario 
when E = E*. As seen from Fig. (|22p . the field sits still on top of the hill while the s field 
oscillates in a finite range. The maximum size of the universe is a finite number determined 
by the constant value of the <fi field. 



D. The case of c < 



There are also geodesically complete solutions when c < and b > which we will outline 
very briefly. The 0, s Hamiltonians are the same as before, as in Eq. ( 1601162]) . with K > 0. 
There exist two interesting classes of geodesically complete solutions that satisfy 2 — s 2 > 0, 
which can occur when c < and b > 0. This happens when the corresponding potentials 
V ((f)) , V (s) take the form in Figs. (25,26). 




The first case is depicted in Fig. (25) when < E < V max (s) = -^rr, and V max (0) = ^ is 
higher, i.e. with b < \c\ . Then s oscillates in the region of a false vacuum, while <fi oscillates 
outside of the false vacuum all the way to infinity. This is similar to the finite bounce 
solution we discussed in section f ]IV CP and figures (I23ll24p . For completeness we record the 
solution 

/ KT 2 _ I /_ x m c (E) = (KT?(E)-l) 
s - (r) = V _f_ sn — \m c , / , \ -V* , (82) 
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and 



I—, out 
7,b>0 



y/KT t 



2 

out 



2bT, 



out an 



T+TQ 
Tout 



"^OUt 



m out (E) = -\{KT 2 out (E)-l) 
T out (E) = [K 2 - WbE) 



-1/4 



(83) 



As in the previous case of the bounce, this solution also represents cyclic bounces at finite 
minimum sizes of the universe, similar to those in Figs.( l23"f24"]) . It occurs for all values of 
b > and c < provided b < |c| , and provided K is large enough so that the curvature terms 
Kef) 2 , Ks 2 are able to compete with the potential term 60 4 , as 4 . The parameters (c, b, K) do 
not need to satisfy any quantization conditions. Also, the synchronization of the relative 
phase is not needed, hence we have allowed the additional integration constant To in the 
solution in Eq.(l83l). 

The second case is depicted in Fig. (26) when < E < V max (s) , where V mSLX (s) = 
and Knax (0) = with b > |c| . Then s oscillates in the region of a false vacuum, while 
the behavior of depends on whether the energy is low (below V max (0)) or high (between 
Knax (0) and Vm ax (s)). The analytic solutions are similar to the ones discussed above, 
except that now c < 0, and the energy E is limited to the region < E < V ma _ x (s) . For 
completeness we list the geodesically complete solutions that satisfy <fi 2 — s 2 > 0, together 
with the quantization condition for their periods 



high 



S 7 ,c<0 l r J 



v , , — sn I 7fr\m c 



m c (E) = (KT 2 (E) - 1) 



8feT + l+cn(^|m H 



</>7,fe>0 ( T ) = 

T + Q (m+) = 2nT c Q (m, 



T c (E) : 

m+ (E) 
T + (E) = 



2 ~ 



¥-4|c|E 



-1/2 



-,\ + KT 2 + {E) 
{<o4bEY l/4 



(84) 



and 



low 



7,c< 



l,b>0 



o(r) 



y/KT^-l 



sn 



(r) = y^sn 



k T b Q (m b ) = nT c Q (m c ) 



m, 



"i'6 



\m b 



m c (E) = (KT 2 (E) — 1 

^ (#) = 
m 6 (S) = (A-.'/;f(A') - I 

T b (£) = 



¥-4|c|£ 



-1/2 



-1/2 



(85) 
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In the low energy level it is also possible for to oscillate on the outside of the false vacuum. 
In that case the solution is given by 



7,c<0 



(r) = v ,_4_ sn [ 7^\m, 



2\c\T c 



T ' \ c I ' 



m c (E) = (KT^(E)-l) 



low : < 



l+,OUt 

7,6>0 



T 



T c (E) = 
m out (£) = -\ {KTl t 



-1/2 



(86) 



* ^(^l" 10 "*) T ou4 (£) ee (K 2 - 166£7)' 



-1/4 



In this case there is no quantization condition on the parameters, and there is an additional 
integration constant r which is arbitrary. 

V. THE OPEN (k = -1) FRW UNIVERSE 



For k = — 1, the dynamics of 0, s is described by the Hamiltonians in Eqs. 
K < 0. The corresponding potentials V (0) , V (s) are plotted in Fig. (127)) . 
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FIG. 27: The open FRW universe, k < 0. 



For the negative energy level E s = E^ < 0, there cannot exist geodesically complete 
solutions that satisfy 2 (r) — s 2 (r) > at all times because of the following argument. The 
energy E s = E^ must be above the minimum of V (s) at the double well to have a solution 
for s. Then the quantity 2 — s 2 always changes sign, whether b > or b < 0. For example, 
suppose s is at the minimum of V (s) , so the solution s (r) is just a constant. But (r) 
oscillates between two turning points m i n (E) < (r) < max (E) , while sometimes 2 — s 2 
is positive and sometimes it is negative. The figure above is drawn for the case |6| < c. If 
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we take |6| > c then V (s) will be like the dashed curve and while V (0) will be like the solid 
thin curve (for b < 0) or the solid thick curve (for b > 0). In these cases again there are no 
solutions such that 2 — s 2 remains positive at all times. 

For the positive energy level E s = > 0, there are geodesically complete solutions that 
satisfy <fi 2 — s 2 > at all times. In fact this case is formally identical to the case discussed 
in section (1IV Aj) . In the present case K < 0, c > 0, while b can have either sign. The exact 
solutions are parallel to Eqs. (164(166]) except for replacing K — — \K\ . Thus, the solutions are 



l l-K^f m {i\ m ') m s {E) = \{l+\K\T 2 {E)) 
MTj V 8cT ° dnfelm.)' T s (E) = (l6cE + K 2 )-^ 1 J 

while cf)^ (r) has the following expressions for b > 

v - m + (E) = \ - \K\T\[E) 
, 7)6 >o (r) = X \ WlTTZJT^ V T+ (E) = {tebEY^ ' (88) 



or b < 



1 sn \%\ 




1 


8bT ll + cn(^\ 





, X 2 T 4 g n(^|m_) m_ (E) = |(l + \K\ T 2 (E)) 

V 8 ' 6 I T2 dn(^y T_(E) = (16 \b\E + Kr 1/4 

As they stand these solutions do not yet satisfy the requirement <fi 2 (r) > s 2 (r) at all 
times. This can be satisfied only by requiring the period of to be a multiple integer of the 
period of s. The analytic expression for this conditions is, as before 

b > : T+Q (m+) = 2nT s Q (m s ) , n = 1, 2, 3, • • • (90) 
b < : T_Q (m_) = nT s Q (m s ) , n = 1, 2, 3, • • • (91) 

Note that in the limit K — > these solutions reduce to the solutions for the flat case with 
b > and 6 < 0. We will not discuss them in any more detail here since this K < case is 
similar to the previous discussion for both K > and K = 0. 

VI. SUMMARY AND OUTLOOK 

We have thoroughly analyzed a simple model of a scalar field interacting with gravity in 
3+1 dimensions. The model was derived from 2T-gravity in 4+2 dimensions as the "3+1 
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dimensional conformal shadow" [4|[5| [6| and can also be constructed^ in the colliding branes 
scenario [2] in 4+1 dimensions using the worldbrane notions [7] inspired by M-theory 

An essential feature of the model in 3+1 dimensions is an underlying local conformal 
symmetry (Weyl symmetry) exhibited in the action of Eq.(JS]). There is no fundamental 
gravitational constant in this model, but instead there is a gauge dependent dynamical 
"gravitational parameter" which plays precisely the role of the gravitational constant when 
the Weyl symmetry is gauge fixed to the Einstein frame, thus agreeing with the standard 
form of Einstein's gravity and its interactions with matter. 

This raised the question of whether the dynamics could force the "gravitational param- 
eter" ((f) 2 (a+) — s 2 (a+)) 1 to change sign in some patches of space-time where antigravity 
would emerge, and whether the existence of such patches could have observational conse- 
quences in our current universe in the context of cosmology or otherwise. This is the question 
that motivated our investigation. 

We found out that generically the gravitational parameter does change sign dynamically, 
and that this change of sign is not a gauge artifact since the gauge invariant quantity 
(1 — s 2 /4> 2 ) can be used to monitor the sign change. So our model indicates that patches of 
antigravity could exist, but we have also found that such patches cannot be reached from 
our current universe without going through singularities, such as the big bang or big crunch 
(perhaps others, such as black holes as well). 

We have not yet answered what the physical implications of this phenomenon may be 
for our current universe, but instead we have limited the current investigation to finding 
and classifying those classical solutions in the context of cosmology that are geodesically 
complete for all times. By all times we mean that one must go beyond the gauge dependent 
definition of "time" and instead seek geodesically complete solutions in all possible choices 
of "time" , thus being able to connect information before and after singularities. In this way 
one is not limited to only some patches while declaring ignorance about other patches 12 . 

We found that the conformal time r is a good evolution parameter for our purpose, so 



12 Of course, we must expect modifications of the classical equations due to quantum gravity especially near 
singularities. We may need to wait a long time before one is sure about what quantum gravity really 
is. Given this cloudiness of our knowledge at the present time, we feel that the geodesically complete 
approach we are pursuing here will still be relevant, and perhaps even provide guidance to clarify the 
issues in future research. 
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we analyzed the solutions for all values of r from minus infinity to plus infinity. In this 
way we learned that the gauge invariant quantity (1 — s 2 (r) /0 2 (r)) oscillates back and 
forth between patches where it is positive and negative (namely gravity /antigravity), and 
this information is carried smoothly through the singularities. In fact we learned that the 
point in time where there is a singularity in the Einstein frame (divergent scalar curvature) 
does not look like a singularity at all in other convenient gauge choices of the Weyl gauge 
symmetry. 

This paper was focussed on finding and classifying the complete subset of all classical 
cosmological solutions for which (1 — s 2 (r) /<p 2 (r)) never changes sign for all times. Thus, 
for the classical solutions exhibited in this paper the universe remains always in the gravity 
patch, never shifting to antigravity. The universe starts expanding with a big bang, but 
eventually it turns around (at a finite or infinite size, depending on the sign of the parameter 
b) and begins to contract, leading to a big crunch. But this is followed with the same periodic 
pattern of a big bang, turnaround, big crunch, again and again indefinitely to the future as 
well as to the past. 

There are such cyclic solutions in which the universe never contracts to zero size, and 
bounces back after contracting to a finite size. So in these finite bounce solutions the universe 
never hits a curvature singularity (as defined in the Einstein frame). These solutions are 
possible for the closed universe. The finite bounce solutions exist for a large range of the 
parameters b, c, K that define the model, but one of the integration constants E (which 
amounts to energy initial conditions for the scalar fields s or 0) must lie in a certain range 
defined by the parameters b, c, K, while the other integration constant tq is arbitrary. 

There are also cyclic solutions in which the universe contracts to zero size periodically, 
thus hitting the curvature singularity (in the Einstein frame) at the big crunches/bangs. 
These cyclic solutions occur for the flat, open and closed universes, but only if some initial 
conditions for the 0, s fields are synchronized (the integration parameter To set to To = 0) 
and a quantization condition is imposed on a combination of the parameters b, c, K and the 
integration parameter E. Thus not every model is capable of yielding geodesically complete 
cyclic solutions, as illustrated clearly in the case of the flat universe. 

Evidently the next stage of this research is to analyze what happens to these solutions 



under perturbations. This is the topic of our next paper in Ref. 



11] where in addition to 



curvature K, radiation and anisotropy are included both at the classical and quantum (in 
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the sense of the Wheeler-deWitt equation) levels. According to our current understanding, 
very similar geodesically complete solutions exist in the presence of these perturbations. The 
question of the physics of antigravity and its effect on our current era of cosmology is an 
interesting topic that we intend to pursue as a natural evolution of the present discussion. 
We hope to report on these details in the near future. 
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